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The scattering of relativistic Dirac particles by a Coulomb field 
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Q^[ tained as a partial wave series. For small Z the series can be summed up 

-t— ' ' 

C ' approximately to give a closed form. The result, though being aproximate, 

qh| exhibites some nonperturbative feature and cannot be obtained from per- 



1 



turbative quantum electrodynamics at the tree level. 



Keywords: relativistic scattering, Coulomb field, two dimensions 
PACS number(s): 03.65.Pm, ll.80.Et 



*published in Phys. Lett. A 260 (1999) 17-23. ©1999 Elsevier Science B.V. 
'fE-mail addresses: qg_lin@163.net, stdp@zsu.edu.cn 
■t Mailing address 

1 



The scattering of particles by a Coulomb field is an important subject in quantum mechanics 
as well as in classical mechanics, because the Coulomb field describes real physical interations 
between charged particles, say, electrons and nucleuses. In the nonrelativistic case, exact 
results of the scattering amplitude or the differential cross section in closed forms can be easily 
obtained by solving the Schrodinger equation in parabolic coordinates [1, 2]. This has been well 
known for a long time. For fast electrons (or positrons) scattered by a nucleus, however, the 
relativistic effect becomes significant and the Dirac equation has to be employed. In this case 
the problem is rather complicated and exact results can be obtained only in open forms, that is, 
in partial wave series. For light nucleuses, the series can be summed up approximately to give a 
closed form for the differential cross sectiion [3, 4]. This includes the nonrelativistic result and 
a relativistic correction. It is the same as the lowest-order result of quantum electrodynamics 
(QED) [5]. Higher-order corrections involve more difficulty. The problem has been investigated 
by many authors over several decades [3-10]. (More references on the subject can be found in 
Refs. [4] and [11].) 

Mathematically one can consider the Coulomb scattering problem in general space dimen- 
sion d > 2. It has been shown that the (nonrelativistic) quantum-mechanical differential cross 
section coincides with the classical one in li = 3 only [12]. Among the various space dimensions 
other than three, the two-dimensional case seems physically relevant. In the recent years, there 
has been wide interest in lower-dimensional field theories and condensed matter physics. The 
Coulomb field may be of interest in some two-dimensional physical models [13, 14] which may 
describe the interaction of Chern-Simons vortex solitons [15-22]. The nonrelativistic scatter- 
ing by a Coulomb field in two dimensions has been solved exactly [23]. The result exhibits 
new features different from those in three dimensions. For example, the quantum-mechanical 
differential cross section is different from the classical one as mentioned above, and the exact 
result is different from the Born approximation. Thus the problem is of some interest in itself 
as well. The purpose of this note is to extend the nonrelativistic result to the relativistic case. 

Consider the scattering of a fast electron (position) with mass fi and electric charge — e 
(e) by a nucleus with charge Ze. When the incident velocity is comparable with c, the speed 
of light in vacuum, the process should be described by the relativistic Dirac equation. We 
begin with the stationary Dirac equation of the electron (position) in the Coulomb field of the 
nucleus: 

Hip = Eij, (la) 

where H is the Hamiltonian 

H = CO. ■ p + 7°Aic^ , (lb) 

r 

where the last term is the Coulomb potential between the electron (positron) and the nucleus, 
K, = ±Ze^ where the plus (minus) sign corresponds to an electron (positron), p is the mo- 
mentum of the incident electron (positron), a = 7^7, and 7^* = (7°, 7) (/i = 0, 1, 2) are Dirac 
matrices satisfying {7^,7^^} = 2g^'^ where g'^'^ = diag(l, —1, —1). In two dimensions the Dirac 
matrices can be realized by 2 x 2 matrices: 

7° = a^ ^' = ^a\ f = ^a^ (2) 

where the a's are Pauli matrices. Thus ip is a two-component spinorial wave function. In two 
dimensions the orbital angular momentum has only one component 

L = e'^xV = -th^, (3) 



where e*-' is antisymmetric in its indices and e^^ = 1, ^ and r used above constitute the polar 
coordinates (r, 6) on the xy plane. It is easy to show that [L, H] ^ even for a free particle. 
However, one can define another operator 



5=^e^V7^ 



(4) 



and show that 



J = L + S 



(5) 



is a conserved operator, that is, [J, H] = 0. Thus S may be regarded as the spin operator and 
J the total angular momentum. J is conserved in any central field. 

The Dirac equation (1) can be solved in the polar coordinates by separation of variables. 
Bound-state solutions (for k, = Ze^) have been studied in detail [24,25]. We are now considering 
scattering solutions with E > fj,c^. We use the representation (2) in which S = ha^/2. Let 



ij,{r,9) 
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(6) 



It is easy to show that Jipj = jhipj. Thus j is a good quantum number. Substituting this 
expression into Eq. (1) we obtain two coupled ordinary differential equations for the radial 
wave functions: 

f--f + hg + ^g = 0, 
ar r r 



where 



— + -g-k2f--f = 0, 
ar r r 



E + ac^ , E - ac^ k 

^1 = ^ ' ^2 = T , 1 = T- 

nc he he 



he 



(7a) 
(7b) 



Then we introduce the new variable 



p = kr, k = \Jk1k2 



^E^ - ^2c4 

he 



and two new functions u{p), v{p) through the definition 

m = le^Xp) + v{p)], g{r) = -'^^/Xp) " Hp)], 



to recast Eq. (7) into the form 



(9) 



(10) 



du i[3 "iP + J 

u V = 0, 

dp p p 

dv id id' — j 

— + 2tv + —v + — -u = 0, 

dp p p 



flla) 



(lib) 



where 






2 VA:^ V^i "^V c- 



P'-i{\l^-r^\-N^-% (12b) 



where Vc is the classical velocity of the incident particle. This is more convenient. Indeed, one 
can eliminate v immediately to obtain an equation for u alone: 

p^ + (l + 2>p)^+f2,5-ll^). = 0. (13) 

ap^ dp \ P J 

where we have used /3^ — /3'^ = 7^. From Eq. (8), 7 ~ ±Z/137. If Z is not very large, say, 
Z < 68, we have I7I < |. Then for any j, the solution of Eq. (13) is well behaved at the 
origin. Let 

u{p) = p'w{p), s = ^f - 72. (14) 

Then we have for w the equation 

p'L^ + {2s + 1 + 2ip)^ + 2{P + is)w = 0. (15) 

ap"^ dp 

This is familiar. The solution that is well behaved at the origin is w{p) = ^{s—if3, 2s+l, —2ip), 
where <l>(a, b, z) is the confluent hypergeometric function. So we have 

Uj{p) = ajp''^{s - if3, 2s + 1, -2ip), (16a) 

where aj is a constant, and the subscript j of Uj that is omitted above has been recovered. 
Substituting this solution into Eq. (11a) we have 

^Ap) = (^.^T^P'Hs -tf3 + l,2s + 1, -2tp), (16b) 

where we have used the formula 

z— + a $(a, b, z) = a^{a + 1, b, z), 

which can be obtained from other relations given in mathematical handbooks [26]. It should 
be remarked that ipj is slightly singular at the origin when j = ±^. However, the integral of 
i^ji^j over any finite volumn converges and the solution is acceptable. We take 



a, = A2-^ij + tf3')^^^^-J^ exp f ^ + imr, - i— + i-\ , (17a) 



where m = j — ^ and 



H^v. 



then we have for the radial wave functions the asymptotic forms 
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fj{r) — * Ai"^ exp{irij) cos ikr + f3ln 2kr 7 + ^i 



r — > 00, 



gj{r) -^ \ —Ai"^ exp{ir]j) sin ( kr + p\n2kr ^ Vj j ^ ''" ^ ^^? 

to the lowest order, where rij are phase shifts defined through 

, T(s — i[3) 
exp{2tr]j) = {j + 1(3 ) j^, . exp(^j7r - isti). 



(18a) 
(18b) 

(19) 



The asymptotic form for ip = J^j'^jj where the summation is taken over all j, turns out to be 



V^ ^ V'in + V^sc, r ^ 00, 



where 



^ii 



l^E- 



lE-fic' 



^i-, 



\\ 2E J 



with 



and 



V^i; 



+00 
m=—oo 



-^ COS Ikr + /5 In 2A;r - — - -') e'^' 
TTkr \ ' 2 4/ 



/ . If + uc" \ 



A 



■ exp{ikr + i(3 In 2kr)f{9) 



2E 



IE - iic^ 



2E I 



with 



(20) 



(21a) 



(21b) 



(22) 



/(^) = \/^Ee^P(^%)sin%e^"' 



^pl^ik 



^[exp(2i77j) - l]e 



imd 



(23) 



In order that ip represents a correct scattering solution, one should show that ipi^ represents an 
incident wave and ipse represents a scattered one. First, it is easy to realize that (^jn is equal to 
exp[i{kr + (3ln2kr) cos6] (for very large r). Thus ^p*n^Phi = 1- The probability current density 
associated with a solution ip is 



j = ctjj'^atjj. 
For ipi^, it is easy to show that 

Jin ^c^xj 

where e^, is the unit vector in the x direction. For ipse we have 



(24) 



(25) 
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jsc = ^-^^v.er, (26) 

r 

where e^ is the unit vector in the radial direction. Thus ip is indeed a correct scattering 
solution, f{6) is the scattering amplitude, and the differential cross section is given by 



a 



{9) = \f{9)\'. (27) 



The choice of aj in Eq. (17) is thereby proved to be appropriate. The results (23) and (27) have 
the same forms as those for a short-range central field [27], but here the asymptotic forms (18) 
where the phase shifts occur involve the logarithmic distortion which is typical for a Coulomb 
field. Since Ej e'""^ = 27r6{e), we have for ^ 7^ 

/(^) = --l=^exp(2^r^,)e^"'- (28) 



V27Tk 



J 



This result with exp(2ir7j) given by Eq. (19) is exact. It is similar to but simpler than the 
corresponding result in three dimensions. However, it is very difficult to sum up the above 
partial wave series in the general case. Closed results can be obtained approximately only for 
light nucleus (small Z). 

Consider light nucleuses which may be more interesting in pratice. With 7 ^ Z/137 and, 
say, Z < 5, we have 7^ ^ 1. Then we may approximately replace s by \j\ in Eq. (19). Note 
that 13 also depends on 7 and we do not make approximation with it. Thus the result will 
possess some nonperturbative features. We have 

, Vim + 112 -id) ,^ ^,J{m+\l2-iS) , , , , 

^^ '^' T{m + 1/2 + il3) ^' ' 'T{m + Z/2 + if3y ^-^ ^' ^ ^ 

, V{\m\ + 1/2 -i(3) ,^ ^,J(\m\-\/2-i^) , , , ,, 

It can be shown that the first term in either equation equals exp(2i(5|m|) where b\m\ is the 
nonrelativistic phase shift of the mth partial wave. The second term is a relativistic correction 
which vanishes when v^/c -^ 0. Substituting Eq. (29) into Eq. (28) we have 

/(^) = /o(^) + /i(^), (30) 

where 

^Y{\/2-i/5) 



Me) 



V2^ 



-F{l,l/2-i(3,l/2 + i(3,e'^ 



ni/2+t(3y 

(31) 



^W^r''^'''^'-''^'^'^''^''^'\ 



is the nonrelativistic partial wave result and 

M9) = -^^I|1Z|^[F(1, 1/2 - tP, 3/2 + tp, e'') - e-^^F(l, 1/2 - tP, 3/2 + tP, e"^^)] 

(32) 



is the relativistic correction. In these equations -F(ai, 02, &i, z) is the hypergeometric function. 
Using their functional relations (cf Eqs. (e.6) and (e.4) in Ref. [2]), it can be shown that 



W) = -I 



r(l/2 - i(3) exp(i/3 In sin^ 6/2) 



nm 



2ks\ne/2 



(33) 



and 



h{0) 



r(i/2-^/3) 

T{i(3) 



^-ie/2 



exp(i/51nsin2^/2) 



2k 



(34) 



where < ^ < 27r. It should be remarked that both terms in Eq. (31) diverge for all 9 [26]. The 
final result (33) is well defined everywhere except at ^ = 0, however. Therefore the calculation 
that leads to Eq. (33) may also be regarded as a regularization procedure. The validity of this 
procedure is obvious because the final result (33) coincides with the nonrelativistic one obtained 
from exact solutions [23]. On the other hand, both terms in Eq. (32) converge everywhere 
except at ^ = [26]. Thus the calculation of fi{0) involves nothing special. Collecting the 
above results we arrive at 



f{d) = -I 



r(l/2 - ip) exp(i/3 In sin^ 9/2) 



2 A; sin ^/2 



1 - le-^'/^ 
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sm 



1 - 






(35) 



We have dropped a 6{9) term in Eq. (28), since it does not contribute when 9 y^ 0. Here we 
have a result more singular than S{9) when ^ — ^ 0. Thus the S{9) term can indeed be dropped 
everywhere and the above expression is enough. It is easy to find that 



^io) = \fm 



/Stanh/Svr 
2k sin^ 9/2 



1- 



9' 



■ sm 



(36) 



Two remarks about the result: First, we have not made any approximation in regard to 
the incident velocity, so the result is valid for high energy collision. It is obvious that the 
relativistic correction becomes significant when Vc is comparable with c. Second, though the 
above result holds for small 7 only, it involves a nonperturbative factor tanh/97r (note that 
/5 oc 7). Thus the result cannot be obtained from perturbative QED at the tree level. This 
is different from the case in three dimensions where the result corresponding to Eq. (36) 
coincides with the lowest-order contribution of perturbative QED. The reason may be that 
in three dimensions the Born approximation (a perturbative result) already gives the correct 
result in the nonrelativistic limit. In two dimensions this is not true. If one expands exp{2ir]j) 
to a higher order in 7 by expanding s and if these higher-order corrections can be summed up 
analytically to give a higher-order correction f2{9) to f{9), then it might be expected that f2{9) 
also possesses the above nonperturbative feature since no approximation is made in regard to 
f3. The higher-order corrections will be studied subsequently. Finally we express a{9) in terms 
of pure classical quantities: 



a{9) 



K,tanh.{7rK,/?iVc 
2fiv^ sin^ 9/2 



1 



-?r sm - 




(360 



Here the first factor is the exact nonrelativistic result. The subsequent ones are due to the 
relativistic effect. They are very similar to those in three dimensions. 
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In conclusion, we have calculated the scattering amplitude for fast electrons (positrons) 
scattered by a nucleus in two dimensions. The exact result is given by Eq. (23) in partial wave 
series. For light nucleuses, the series can be summed up approximately to give a closed result 
(35). The differential cross section is given by Eq. (36) or (36'). Though being approximate, 
the results exhibit some nonperturbative feature and cannot be obtained from the lowest-order 
contribution of perturbative QED. This feature is not manifest in the corresponding result in 
three dimensions. 
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supported by the National Natural Science Foundation of China. 
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